Two-photon decay of heavy hadron molecules 
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We discuss the two-photon decay width of heavy hadron molecules and study the dependence on 
the constituent meson masses and on the binding energy. In addition finite size effects due to the 
extended structure of the bound state are shown to have a strong infiuence on the predictions for 
this decay width. 
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Conventional quark-antiquark constituent quark 
model descriptions of the meson spectrum are certainly 
incomplete and sometimes inadequate to match ex- 
perimental observation. As a complementary or even 
alternative approach new structure assumptions have 
been invoked in the past. This concerns for example 
glueballs, quark-gluon hybrids, tetraquark configurations 
or even an admixture of some of these states vifith an 
ordinary quark-antiquark meson as suggested by direct 
applications of (as in lattice simulations) or in models 
motivated by the underlying theory of QCD [l|, 01 • But 
also conventional configurations can populate the meson 
spectrum as quasi-nuclear nucleon-antinucleon bound 
states j3| or hadronic molecules, bound states of tviro or 
more mesons [4-6]. 

The idea that mesons can also form mesonic bound 
states traces back to the original work of 0, Q with a 
first systematic study for the meson spectrum given in [9| . 
The concept of hadronic molecules is often invoked when 
observed meson masses are rather close to a respective 
two-meson threshold and when dynamical arguments can 
be found which suggest binding in these systems. Exper- 
imental candidates for these type of states are especially 
discussed in light of the experimental advances in the 
heavy flavor sector. A prominent example is the X(3872) 
where a molecular interpretation is suggested lid] . 

A full interpretation of an observed meson resonance 
is not only based on a consistent description of its mass 
and J^'^(/*^) quantum numbers, but also on a possible 
explanation of the observed decay and production modes 
of this respective state. Here the electromagnetic inter- 
action plays a leading role since it is well understood even 
on a effective hadronic level and can be treated perturba- 
tively. In particular, two-photon or, in general, radiative 
decays are considered diagnostic tools which are sensitive 
to the iimer structure of the short-lived resonance @ . 

The purpose of the present manuscript is the discus- 
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sion of the two-photon decays of heavy hadron molecules. 
Our aim is to work out the influence of the binding energy 
and the constituent meson masses on quantities like cou- 
plings (of the molecular state to its constituent mesons) 
and electromagnetic decay amplitudes. Under certain 
circumstances, as for a local theory, the results are ana- 
lytical and the dependence on the above mentioned quan- 
tities can be shown explicitly. Especially ratios of two- 
photon decay rates of bound states, also in different fla- 
vor sectors, provide a simple and clear estimate whether 
a molecular structure is likely or not. 

For the physical or nonlocal case, which effectively 
models finite size effects of hadronic molecules, results 
can deviate sizably from the limiting case of a local the- 
ory. Here we find that, in general, the two-photon decay 
properties of heavy bound states are more sensitive to fi- 
nite size effects than hadronic molecules of light mesons. 

Several heavy meson and baryon states are considered 
good candidates for a molecular structure among them 
the L»*o(2317), i:»^i(2460), X(3872), y(3940), r(4140), 
F(4660) and Ac(2940). For a review of the situation in 
the heavy meson sector we refer to [l|, ]% UM ■ 

For the present study of two-photon decays of heavy 
meson molecules we choose the framework of a phe - 
nomenological Lagrangian approach developed in [5|, Il2l | . 
The bound state structure of hadronic molecules is set 
up by the compositeness condition [l^, [l3] which im- 
plies that the renormalization constant Z of the hadron 
wave function is set equal to zero or that the hadron 
exists as a bound state of its constituents. Decay pro- 
cesses of hadronic bound states are described by the 
coupling of the final state particles via meson-loop dia- 
grams to the constituents of the molecular state [a. la. Il4| . 
First calculations for the two-photon decay widths of 
heavy hadron molecules as F(3940) = {D*D*'^} and 
r(4140) = {D;dP} have been performed in [l|. 

In this note we concentrate on hadronic bound states in 
the meson sector. We further restrict to S-wa.ve hadron 
molecules with quantum numbers J^^ = 0^^ whose 
constituents are pseudoscalar or scalar charmed/bottom 
mesons. The coupling of the scalar molecular state H 
to its constituents is expressed by the phenomenological 



Lagrangian: 



Cnix) = gHH{x)JH{x) 
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where gn is the couphng constant, H is the field repre- 
senting the hadronic molecule and Jh{x) = (/)^(a;)(/)(x) 
is the current carrying the quantum numbers of the 
molecule composed of the constituent mesons identified 
with the field (j). A nonlocal extension, which corresponds 
to the physical case of a hadronic molecule with a finite 
size, can easily be achieved by insertion of a vertex func- 
tion $/f (y^) into the current Jh{x): 

Jh{x)^ fd*y<^H{y^)^H^ + y/2)H^~y/2) (2) 



where y is the relative Jacobi coordinate. Here ^niv^) 
is the correlation function describing the distribution of 
the constituents inside the molecular state H. In the 
present evaluations we use a Gaussian form with $//(y^): 
^h{pI:/J^hJ = exp(-p|/A|^.), where pe is the Eu- 
clidean Jacobi momentum and A//, is a size parameter. 
The coupling constan t qh is determined by the com- 
positeness condition [il,[i3, where Z — l—gj^Y.'{AP) — 
0. Here, ^'(M^) = dT,(p'^)/dp'^\p2^M^ is the derivative 
of the mass operator S(p^) illustrated in Fig. [H where 
M is the mass of the hadronic molecule. The generic ex- 
pression for the mass operator of a bound state of two 
pseudoscalars is e.g. given in [6], 
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where S{k) = {m? — k"^ ~it)~^ denotes the free propagator 
of the constituent meson. 

In the following we first restrict to the case of a local 
coupling of the H state to its respective constituents. In 
case of a local interaction the vertex function is replaced 
by $//(— fc^) = exp(fc^/A|^.) = 1. Then the coupling gn 
can be expressed in the analytical form [g, [l^| : 
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where C, — M/(2m), /?(C) — arcsin(C)/C and m is the 
mass of the constituent meson. 

Inclusion of the electromagnetic interaction in a gauge 
invariant way is discussed in our previous works [y, [la] ■ 
In the local approximation we deal with the two diagrams 
of Fig. [21 The extension to the nonlocal case leads to 
further diagrams, which need to be included in order to 
guarantee gauge invariance (a full discussion can be found 
ini[i6]). 

The gauge invariant form of the matrix element for the 
radiative transition in case of real photons reads as: 



M,, 
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where qi and 92 are the 4-momenta of the photons. For 
the local case the expression for the effective coupling 




FIG. 1: Mass operator of the hadronic molecule H = {(j>4>^) 
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FIG. 2: Diagrams contributing to the H — > 77 decay. 



constant gn-yy is given by the simple form: 
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However, this expression becomes more complicated and 
can only be solved numerically when including finite size 
effects (see e.g. Appendix of 6|). The radiative decay 



width of the 0^ 
formula 



molecule is calculated according to the 



T[H ^ 77) 
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In the local limit we can write the decay width in terms 
of the quantity C, by using Eq. ^ 



a 



with 



F(ff^77) = — m/(C), 

ZTT 
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Since in the case of hadronic bound states the binding 
energy e is small in comparison to the masses of the con- 
stituent mesons m we can perform an expansion of gu, 
gwy-i and T{H — >■ 77) in x = e/(2?7i), where Q = 1 — x. 
In order to guarantee an accurate approximation of the 
observables we need to include the leading (LO) and 
next-to-leading order (NLO) terms in the expansion of 
gn and gnj-y Therefore, the expansion of T{H -^ 77) 
includes three terms up to the next-to-next-to-leading 
order (NNLO) contribution. 

The corresponding x-expansions of the quantities of 
interest are given by: 
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In the next step we consider bound states of the pseu- 
doscalar D, Dg and B mesons which have the following 
structure 
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Note, the existence of bound states of two heavy pseu- 
doscalar mesons was proposed before in Ref. [171] where 
hadronic molecules are dynamically generated in a cou- 
pled channel formalism. Based on the identifications of 
Eq. ([T2l) additional flavor factors have to be considered 
in Eq. (flUl) which leads to: 
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where Xi = e/(2mi), CHij-y = l/\/2 for i = D,B and 1 for 
i = Ds- From last equation follows the ratio of the two- 
photon widths of two different molecular states which is 
characterized by the respective constituent masses and 
binding energies. At leading order we deal with the sim- 
ple expression: 



TJHb -^ 77) ^ / mBeg y/^ 
r(i/D ^ 77) \mDeD-' 



(14) 



Varying the binding energy e from 10 to 100 MeV which 
(besides the exception of the X(3872)) are typical values 
in the heavy meson sector, the two-photon decay widths 
are evaluated as 



r(iJz5^77) = 0.25- 1.19 keV, 

r(7JD, ^77) = 1.31- 2.58 keV, 

r(ffB^77) = 1.19- 2.78 keV. 



(15) 



Since T{Hi — > 77) oc y/me the smaller value of the bind- 
ing energy also results in a smaller radiative decay width. 
Up to now the results for the coupling and the two- 
photon decay widths were given in the local limit. For the 
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FIG. 3: Couplings gHj^-it and QHi-i-f with i — D,Da,B in 
dependence on Ah, for ex = 7.35 MeV and e^ = 10 MeV. 



physically appropriate description of an extended object, 
the hadronic molecule, finite size effects in terms of the 
vertex function have to be included. This will lead to a 
suppression of the couplings and the two-photon decay 
widths, hence the local case presents an upper limit for 
these quantities. In the following we study in addition to 
the influence of mass and binding energy the dependence 
on the size parameter A//. , which models the nonlocality. 

The dependence of the couplings gHi-y-y on the size pa- 
rameter Afji , typically chosen in the range from 1 to 3 
GeV, is illustrated in the logarithmic plots of Figs. [3]- 
m for different values of e ranging from 10 to 100 MeV. 
For comparison we also include the coupling guKi-y of 
/o(980)/ao(980) to two photons with the well determined 
binding energy of e = 7.35 MeV. The local case, i.e. 
A//. — >■ 00, characterizes the asymptotics of the curves. 
The convergence of the coupling gHi-y-y towards the local 
or asymptotic value depends on the constituent meson 
masses. The coupling gHK-y-y i^ almost stable with respect 
to variations of Ah. near 1 GeV. The couplings of the 
heavy hadron molecules are more sensitive to finite size 
effects for values near A/f. w 1 GeV, note that in Figs. [Hi- 
ll] the dependence on A//, is displayed on a log scale. In 
fact, the coupling gnKii reaches 90% of the asymptotic 
value (local approximation) already at A/f . « 0.6 GeV. In 
contrast, the couplings of heavy bound states approach 
90% of the local value at around 4 GeV in case of the D 
and Z3s-meson bound states and at about 9 GeV for the 
even heavier hidden-bottom molecule. 

The behavior of the couplings gHi-y-y, which enter 
quadratically in the expression for the radiative decay 
width, shows that the suppression due to the size param- 
eter A/f . is larger in case of heavier constituents. Phe- 
nomenologically, the values for A//, tend to increase in 
case of heavy bound states. This behavior was already 
observed in earlier analyses in the framework of meson 
molecules [3, 0, [ij b^t also in case of baryons [l3| • But 
still, for reasonable values of Ajy^ heavier systems are 
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FIG. 4: Couplings gHj^-y-y and gu^ 



with i 



D,Ds,B in 



dependence on Ah; for e^ = 7.35 MeV and e^ = 100 MeV. 



From previous discussions it should be clear that finite- 
size effects are quite important for a quantitative deter- 
mination of the radiative decay width of heavy hadron 
molecules. The same observation holds for molecular 
states composed of vector mesons [U] , where we showed 
that for Ahjj,- = ^Ho* =1^2 GeV the radiative widths 
of the molecules y(3940) = {D*D*^ and y(4140) = 
{D;Dp} are of order of 1 keV for A^, = 2 GeV. 

For completeness, we compare the two-photon decays 
of the heavy systems to the one of the light scalars 
/o/oq. The stability of the coupling in case of a KK 
bound state already implies that the radiative decay 
width r(/(980) — > 77) is not sensitive to the cutoff or 
finite-size effects provided A//, is above 0.5 GeV. Here 
we completely agree with the conclusions of Ref. [15|'). 
In particular, in order to reproduce the current data on 
strong and radiative decays of /(980) we fixed the cutoff 
parameter A^^ — ^/ = 1 GeV [y] with 



strongly influenced by finite size effects as reflected in 
the two-photon coupling gHij-y- 

In practice Ah^ should also depend in average on the 
binding energy. For instance, a small binding energy 
would lead to a loose bound state with a more extended 
structure than a strongly bound compact state. There- 
fore A//, should decrease with smaller binding energies, 
increasing the deviation from the local limit even more. 

To quantify the inclusion of finite size effects for the 
two-photon decay widths we also give results for two val- 
ues of Ah,. For A^. = 1 GeV and e = 10 - 100 MeV we 
have 



r(/(980) -^ 77) = 0.25 keV 



(18) 



T{Hd -^ 11) = 0.05 - 0.34 keV , 

T{Hd, -^ 77) = 0.32 - 0.68 keV , 

T{Hb -^ 77) = 0.05 - 0.12 keV . 



(16) 



When increasing A/f . to 2 GeV the results are by a factor 
2—3 larger 



V{Hd^11) = 0.13- 0.73 keV, 

T{Hd, -^ 77) = 0.71 - 1.52 keV , 

T{Hb-^11) = 0.18- 0.44 keV. 



(17) 



These results should be compared to the local limits of 
Eq. (HSl). 



a value which is very close to the result of the local ap- 
proximation: 



r(/(980)^77) =0.29keV. 



(19) 



Similar results are obtained by the molecular approach 
in [15] and |19j. However, for smaller values of Aj < 
0.6 GeV the radiative decay width decreases which means 
that such small values of the cutoff parameter Aj are 
unlikely according to present data. 
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